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Abstract. A theorem, giving necessary and sufficient condition for naked singularity 
formation in spherically symmetric non static spacetimes under hypotheses of physical 
acceptability, is formulated and proved. The theorem relates existence of singular null 
geodesies to existence of regular curves which are super-solutions of the radial null geo- 
desic equation, and allows us to treat all the known examples of naked singularities from a 
unified viewpoint. New examples are also found using this approach, and perspectives are 
discussed. 



1. Introduction 

The validity of the so called Cosmic Censorship Conjecture, i.e. the idea that physically 
acceptable collapsing systems should always form blackholes is, as is well known, seri- 
ously put in doubt by various counterexamples, that is, exact solutions of the Einstein field 
equations exhibiting naked singularities. 

Recently, we re-considered this problem in an effort of studying the validity of cen- 
sorship, at least in spherical symmetry, with the application to the geodesic equations of 
techniques of non-linear o.d.e. without resorting to exact solutions of Einstein's. It is, in 
fact, obvious that exact solutions, despite their obvious interest, cannot give the complete 
answer to this problem due to the non-linearity of the field equations. 

Our first results, presented recently adopted the 'censorist' point of view, in that we 
provided a sufficient condition for black hole formation in spherical symmetry. However, 
the existence of such a condition does not necessarily mean that the room available for the 
Censor is very large (that is, that the hypotheses of the theorem are likely to be realized in 
generic situations). In order to face this problem, we take here the opposite viewpoint and 
obtain a sufficient condition for naked singularity formation (the condition, as we shall see, 
is also trivially shown to be necessary). The theorem allows us to give a simple, complete 
classification of (virtually) all the widespread zoo of naked singularities existing in the 
literature, and to produce new examples as well. 

2. The existence of naked singularities 

We consider a spherically symmetric collapsing object (the matter source can be any 
model compatible with the weak energy condition). The general, spherically symmetric, 
non-static line element in comoving coordinates t, r, 9, ip can be written in terms of three 
functions \, Rof r and t only as follows: 

ds^ = -e^^df + e^^dr^ + R^{de^ + sin^ Od'p^) , 

A fundamental quantity is the mass function m(r, t) defined in such a way that the 
equation R = 2rn spans the boundary of the trapped region, i.e. the region in which 
outgoing null rays re-converge: 

(2.1) m(r, t) = (i?/2) (1 - gf^^d^Rd.R) . 

The curve th{r) defined via R{r, th{r)) ~ 2m(r, th{r)) is called apparent horizon. 

We consider here only those matter configurations, which admit a regular center, and 
we always suppose that the collapse starts from regular initial data on a Cauchy surface 
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(t = 0, say), so that the singularities forming are a genuine outcome of the dynamics. 
If the solution is initially regular, the spacetime can become singular whenever R = 
(focusing singularities) or i?' = (crossing singularities). We consider here only focusing 
singularities. 

The locus of the zeroes of the function R{r,t) defines the singularity curve ts{r) by 
the relation R{r, ts{r)) = 0. Physically, ts{r) is the comoving time at which the shell of 
matter labeled by r becomes singular The singularity forming at r = 0, < = ts(0) is called 
central as opposed to those occurring at r = rp > 0,t = ts{rQ). A singularity cannot 
be naked if it occurs after the formation of the apparent horizon. Since R vanishes at a 
singularity, any naked singularity in spherical symmetry must be massless; at a massless 
singularity the horizon and the singularity form simultaneously. Since regularity of the 
center up to singularity formation requires m(0, t) = OVt < ts{0), the center is always a 
candidate for nakedness. Usually, all other points ("non-central points") of the singularity 
curve lie after the formation of the apparent horizon, i.e. th{r) < ts{r) for r > 0, and the 
singularity is therefore covered. This always happens, for instance, if the radial pressure 
is positive. Then we concentrate here on the central singularity, although our results can 
easily be extended to non-central singularities as well. 

To analyze the causal structure, observe that, if the singularity is visible to nearby ob- 
servers, at least one outgoing null geodesic must exist, that meets the singularity in the 
past. Such a geodesic will be a solution of 



(2.2) 



i(0)-ts(0), 



where (p{r, t) :— e^^'^ . For a problem of this kind, in which the initial point is singular 
(the function Lp is not defined at (0, ts(0))) no general results of existence/non existence 
are known. 

In what follows we ere going to consider sub and super solutions of (2.2). We recall 
that a function j/o(^) is called a subsolution (respectively supersolution) of an ordinary 
differential equation of the kind y' = f{r, y) if it satisfies ?/q < fir, yo) (respectively >). 
In ||3| we have shown the following results: 

2.1. Lemma. If the weak energy condition holds and the radial stress is non-negative then 
the apparent horizon th{r) is a subsolution of (2.2). 

2.2. Theorem. If the weak energy condition holds and dLpjdt < in a neighborhood of 
(0, ts{0)) the singularity is covered. 

2.3. Remark. The above results hold also in presence of radial tensions (i.e. negative radial 
stresses) provided thatp^ satisfies, near the centre, the bound SnR^Pr > — 1. 

We take now the opposite point of view: we search for conditions for naked singularity 
formation. 

2.4. Definition. A curve i+(r) is called a sub-horizon supersolution (SHS) if 
(2.3) t+(0) = th{0), t+{r) < th{r) t'^{r) > (t>{r,t+{r)) Vr > 0. 

2.5. Theorem. The singularity is naked if and only if a SHS exists.. 

Proof. If the singularity is naked the singular geodesic (the solution) is also a supersolu- 
tion. We now proof sufficiency. Take a point (ro,io) in the region S ~ {{r,t) : r > 
0, t+{r) < t < th{r)}. At this point the (regular) Cauchy problem for ip admits a unique 
local solution tg{r). Now the extension of this solution in the past cannot escape from 
S since either it would cross the supersolution from above (leading to a regular geodesic, 
eventually escaping from the central point at a time prior to singularity formation) or it 
would cross the subsolution from below. Thus it must extend back to the singularity with 
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2.6. Remark. It is easy to show that, if a SHS exists, then ^ cannot be strictly negative in 



a neighborhood of (0, to), in accordance with theorem 2.2 



The above analysis is limited to radial null geodesies. However, it can be shown that, 
if no radial null geodesic escapes from the singularity, then no null geodesic escapes at all. 
In other words, we show that a radially censored singularity is censored In the special 
case of dust spacetimes, this result has been was first given by Nolan and Mena []l4|]. 

2.7. Theorem. A radially censored central singularity is censored. 

Proof. Suppose, by contradiction, that the center is radially censored but that there exist a 
singular non radial null geodesies {i{r), say). This curve satisfies to 

dt 



(2.4) — = Je2^-2'' + e-2=^i2/i?2 > ipir, t), 

where is the conserved angular momentum. Thus t{r) is a supersolution of the null 
radial geodesic equation. Obviously, t{r) < th{r) so that t{r) is a SHS. As a consequence 



of Theorem 1.5, there exist a singular null radial geodesic i.e. a contradiction. □ 

3. Old and new examples 

Comoving coordinates r, t are extremely useful in dealing with gravitational collapse 
because of the transparent physical meaning of the comoving time. We are, however, go- 
ing to use in the present section another system of coordinates, the so-called area-radius 



coordinates, which were first introduced by Ori 1 15| to study charged dust, and then suc- 
cessfully applied to other models of gravitational collapse (see e.g. ^^). These coordi- 
nates prove extremely useful for technical purposes. 

In area-radius coordinates the comoving time is replaced by R. The velocity field of the 
material v^^ = £"'^6^ transforms to — e^'^ R5^ and therefore the transformed metric, 
although non diagonal, is still comoving. In a recent paper, we have found the following 
class of solutions of the Einstein field equations in area-radius coordinates [Q: 

/ 2^'\ Y 1 

(3.1) ds^ = - 1 G^dr^ + 2G —dRdr ^dR^ + R^idO^ + sm^0dLp^) . 

\ R J u 

In the above formulae, the two functions '^{r, R) and Y{r, R) are arbitrary (positive) 
functions while 

(3.2) y2 = Y^+'^-^-l. 
and the function G is given in terms of a quadrature: 

1 1 



(3.3) G(r,i?)= / —L_^ilM{r,a)da. 

J p. Y{r,a) or Y[r,r)u(r,r) 

The matter distribution of the source is described in terms of energy density, radial and 
tangential stresses as follows (a comma denotes partial derivative): 



(3.4) 



inR^YuG AttR'^ ' 



(3.5) pr - -7-7^, Pt 



47ri?2 ' ''^ SnRuG \ Y Y^ J Svri? 

Special, very well known sub-cases of the metrics above are: 

1) The dust (Tolman-Bondi) spacetimes. All stresses vanish, energy density equals the 
matter density, and this implies '^^h = and Y — E'^^r- 

2) The general exact solution with vanishing radial stresses. Vanishing of pr implies 
'i' li — 0, while Y depends, in general, also on R. The properties of these solutions have 



been widely discussed in [ 12 1 
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We are now going to show that Theorem all ows for a complete classification of all 
the naked singularities contained in the metrics (3.1). Since in turn such metrics contain all[| 
the special cases for which censorship has already been investigated, we also get a unified 
treatment of all these - quite scattered around - results. In every such cases, the nature 
of the endstates has been obtained after a considerable amount of joint work by several 
authors, sometimes developing original mathematical techniques such as the so called root 



equation technique formalized by Joshi and Dwidedi 1 11 1 . 

The arbitrary functions \1/ and Y must be chosen in such a way that the metric admits a 
regular center prior to singularity formation and that the weak energy condition is satisfied. 
It is not difficult to check that this implies constraints on the first non-vanishing terms of 
the two arbitrary functions near the center It is, however, useful for technical reasons to 
express such constraints in an equivalent way using Ru^ and Y^, as follows: 

(3.6) Ru^ = J2 h^jr'W + h,jr'R' +..., 

i+j=3 i+j=3+p 

(3.7) Y^{r,R) = l + fir' + ^ fc„r*i?^' + . . . . 

i+j=q+l 
j>0 

where ^,30 (that we will call a hereafter) is a strictly positive quantity, / > 2 and dots 
denote higher orders terms. 

It can now be shown that the following Taylor expansion of G{r, 0) near the centre 
follows: 

(3.8) G{r,0) =parP-'^ -br", 
where, defined 



(3.9a) P(T)=^/iyT^ Qir) = - ^ h,, 



i+j=3 i+j=3+p 



(3.9b) S{t)=. n^)= E ^'^'J^'' 

i+j=q+l i+j=3 
j>0 

it is 

Thus, we finally have 

(3.11) G(r,0)=Cr"-i + ..., 

where ^ > and n is a positive integer, namely the smallest between p and q + 1, although 
very special cases - where the two terms balance each other and one must compute higher 
order terms - can conceived. In particular, the second term is related to the acceleration of 
the matter flow lines and vanishes if this does. 
The folUowing can now be proved [B] : 



3.1. Theorem. In the spacetime described by the metric (3.1), a SHS exists - and therefore 
the singularity forming at R = r = is naked - if and only if n — 1, n = 2, or n = 3 and 
S, > a(,c where = 26+i^. 

In what follows, we briefly discuss the cases already known as well as some new ones. 



'As far as we know the unique exceptions are the Vaidia spacetimes |jl||, which however can be easily accom- 
modated in the present formalism, and the case of self-similar perfect fluids, where however the field equations 
become ordinary differential equations and the standard techniques of dynamical systems can be applied Ipp. 
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3.1. Dust clouds. In dust models both vj/ and Y — depend only on r. The acceleration 
vanishes and therefore the index n = p. The solution can be given in explicit form. The 
spectrum of endstates for these models is very well known and was first calculated in full 
generality by Singh and Joshi JlO[|. However, our unifying framework somewhat changes 
the way in which the physical content of this structure must be understood, so that we 
discuss it briefly. 

Introduce the Taylor expansions 



(3.12) 



.3+11 



Considering first the so-called marginally bound case (this case has been recently re- 
analyzed via o.d.e. techniques in Q) Y = 1, we have that the index n = p = h 
and a = 2Fq, P{t) — 2Fo, Q{t) — —2Fn- The central singularity is naked if n = 1, 2, 
censored for n > 4. In the critical case n — 3 the transition is governed by the parameter 
adlTol): 

V7 2F3 2F3 ^ 2Fo^ 

" io 2{2F,f/^ ^ 3(2Fo)3/2 > 3 
In the general (non marginally bound) case, the function Y is no longer constant. As a 
consequence, the expansion of the initial density and velocity mix up to produce the index 
n and the transitional behavior. The singularity is naked if (i^i, /s) 7^ 0, or if (i^i, /s) — 
but {F2t fi) ^ 0, corresponding respectively to the cases n ~ 1,2. The critical case 
occurs when {F^.h) = {F2J4) = 0, but {F3J5) ^ 0. Since a = 2Fo, P(r) = 
2-Fo + /2T, Q{t) — — (2F3 + Jzt), for the singularity to be naked|^it must be 



(2f3 + hr)V^ 
2(2Fo + /2r)3/2 



dr 



A- 

2Fo' 



El 
Fo 



3 is 

2 /2 



h 

2Fo 



F3 

Fo 



h 
h 



^2Jo 

> 3 4c, 



where the function F is defined as 



(3.13) 



ny) 



arcsiiih^— y 



OEv. 



arcsin y'y 

^372 



for y < 0, 
for y = 0, 
for y > 0. 



Finally, if {Fi, /s) = (-^"2, /4) = (-P3, /s) = the singularity is censored. 

To discuss the physical meaning of this structure we start from an interesting example 
given by Dwidedi and Joshi 

Consider a mass distribution which generates a initially homogeneous energy density 
(^'(r) ~ For^). Now, keeping fixed ^, consider different initial velocities. If the cloud 
starts collapsing from rest the resulting spacetime is nothing but the "paradigm" of black- 
hole formation, namely, the Oppenheimer-Snyder solution. However, it suffices a fifth- 
order Y^, say y^(r) = 1 + for^ + f^r^, to change drastically the behavior of the system. 
In fact, the fifth order term generates the critical case. 

What actually happens is, that the mechanism controlling existence of SHS's depends 
only on the value of n, independently from the details of its physical origin. The mecha- 
nism works as a "kibbling machine" looking at just one term of the Taylor expansion of 
just function near the center It accepts any (physically sound) input and returns an output 
which depends on the value of an integer: we can call it a n-machine . 

As we shall see in next example, also the equation of state plays a similar role of "just 
an input" for the n-machine. 



■^In ijl^ a quantity called Q3 plays the role of a, and its expression should be connected - formula (36) in [ [ic| ] 
■ to Q3 = Sav'Tb, with a given above. 
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3.2. Vanishing radial stresses. These solutions are characterized by a mass depending 
only on r, but by a function which contains the contribution of the internal elastic 
energy and thus depends also on R. As a consequence, the metric cannot be given in 
explicit form in general (a notable exception exists however [^). 

For simpUcity, we consider here only the marginally bound case. To describe the data it 
is convenient to consider the Taylor expansion of the mass and of the function g = — 1. 
The constraints of physical acceptability then imply [ 12 1: 



hr'^ir-Rf 
1 + (3kr''{r - Rf 



where fc > 1. Although the solutions are accelerating, it can be easily shown that acceler- 
ation terms do not enter into the first non-vanishing term of formula ( 3.11 ). We thus need 
to consider the Taylor expansion of Ru^, given to lowest orders by 2For^ + 2Fnr^^^ — 
[ikr^ R (r^ + i?^ _ 2rR). Since in general the solution is known only by quadratures, 
when the spectrum of endstates was first calculated in full generality ^ the Taylor expan- 
sion was used under the integral using approximation techniques near the center (delicate 
estimates coming from Lebesgue dominated convergence theorem are required). The end- 
states look like to depend on two parameters n, fc, the second of them coming exclusively 
from the non-trivial equation of state. However, again, if we use the n-machine of Theorem 
3.1 we see that the mechanism has no interest at all in the physical source of the parameter 
n: a naked singularity (rt = 1 , 2) can originate from fc = 1 , 2 or from n = 1,2. The critical 
case always occurs at n = 3 and what changes is only the value of the critical parameter, 
since it must be akn > ^^^c, where^ 



2{2Fo)y^ 



dr 



(2Fo)3/2 



6 61b 



3.3. Acceleration free non-dust solutions. In recent years, spacetimes with cosmologi- 
cal ("lambda") term have attracted a renewed interest both from the astrophysical point of 
view, since recent observations of high-redshift type la supernovae suggest a non-vanishing 
value of lambda, and from the theoretical point of view, after the proposal of the so-called 
Ads-Cft correspondence in string theory. The unique model of gravitational collapse with 
lambda term available so far is the Tolman-Bondi-de Sitter (TBdS) spacetime, describing 
the collapse of spherical dust (within the framework of the present paper, these solutions 
are obtained choosing ^(r, i?) = F{r) + ^R'^)- The spectrum of endstates for these 
solutions has recently been obtained . 

One way to construct new collapsing solutions with lambda term is to choose a function 
of the type ^'(r, R) = F{r) + M (R) R^, with M (R) not necessarily constant as in TBdS. 
It is easy to check that, in order to satisfy the requirements of physical reasonableness, one 
must have *(r, R) = For^ + Fnr^+'' + MqR^ + MkR''^^ + .... Then, if one among 
{Fi , F2 , Ml , M2 ) is non zero the singularity is naked. If otherwise they all vanish, and one 
among (F3, M3) does not, we must compute the parameter a ( 3.10 ): 



(2Fo + 2Mot3)3/2 



dr 



3(2Fo)3/2 



^0 



F3 + M3 



1 



Mo 
Fo 



In nsll a multiplying factor is missing. 
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where F is given in (3.13). The singularity is naked if a > ^^£,c- In all other cases the 
singularity is censored.^ Once again, the physical parameters of the collapse (in this case, 
there is also the cosmological constant) are mixed up by the n-machine. 

4. Discussion 

Starting from the piooneristic work by Eardley and Smarr [ [l6t and Christodoulou [ |l7| 
on dust clouds, the number of papers containing analysis of examples of naked singular- 
ities in spherically symmetric spacetimes can, up to now, be estimated to more than one 
hundred. Authors have tried to make clear the connection between, on one side, the final 
states of the collapse and, on the other side, the choice of the initial data and of the matter 
models. To do this, a mathematical framework was developed (the so called root equation 
technique) which allows investigation of the behavior of solutions of a differential equa- 
tion at a singular point. To be applicable, this framework requires explicit knowledge of 
the solution of the Einstein field equations in which the geodesic motion is studied: it can- 
not be applied if only a part of the field equations are integrated, and can be applied only 
with enormous technical efforts if the solution is known only by quadratures. It is anyway 
obvious that, due to the non-linear character of the Einstein equations one cannot have any 
choice of giving general answer to the censorship problem, even in spherical symmetry, 
using only exact solutions. 

So motivated, we started re-considering the censorship problem in spherical symmetry 
in an effort of going beyond the need for exact solutions. Essentially, the idea was to rely 
as much as possible on the "differential level" - the fact that the metric satisfies to the 
field equations - without resorting to the fortunate case of explicit solutions. This approach 
leads naturally to the application of comparison techniques in o.d.e. since the physical 
properties of the apparent horizon have as mathematical counterpart the fact that this hori- 
zon is a sub-solution of the geodesic motion ||^]. Here we presented a new framework for 
the investigation of censorship, based on such ideas. The framework can be easily applied 
when the exact solution is known only up to quadratures, and this allowed us to give a uni- 
fied view of (virtually) all the examples of naked singularities which are scattered around 
in the literature, as well as to produce new examples. The application of such techniques 
allows, however, also to investigate cases in which the metric is not known at all, such as, 
for instance, the case of barotropic perfect fluids. 
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^In the particular case when all the vanish for > 0, we recover Tolman-Bondi-de Sitter spacetime. 

The condition for nakedness at the transition reduces to > slightly correcting a wrong 

(2Fo) = /2^/^ + l 

value given in 
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